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^STBiNCT 


i:<ifidl vmm and land eov«r infowaatioii Ic very inijortant 
for tho reoemre® planner** H«ne« this study tendsto have 
Imense sodo-aconorale li^xartanee* Thm mm of panchroiaatlc 
ttedlwa m^Xm aerial jhatogpraiiia and 3att<®r the ij»e of 
siaaU, imle pliotograjlia ad3ted ultti sattelit* plctiaras haw 
helped a lot In this sphere* Hi® os® of iahdsat t©s data 
thou^ having a coarse resolMtlon ar© of ws® fear this 
plannliigw i^lng @o m zwMolijtloii data th® mxvacmsf of ^S^% 
has h®eii dhtalned* Certainly idils situation has scop® of 
dswolopnoiit by using th® sensor data that have good msoli** 
tlon as oosparad to the umshT ms* Hi® iis@ of Wmfm 
clai i sl fl s r idth a-prlorl .infcrsatlcsi abcnit th® leisd us® 
pattern cm osrtalaly laprov® the aascorewry of prepared 
Zand m^* Hi® r^p®t4tlv® natumi of ihita anrallsible 
from sattelltes h®]ps in assessing th® tmqporal* spatial 
ehwngas of lend rssouro®® hiiis® Justify tihalr ism* 



CHAPTER 1 


INTRODUCTION 


1*1 Introduction ; 

Lend use and land cover data are essential to 
planners who have to make decisions concerning land 
resource management. Hence they tend to have immense 
economic importance. The term land cover relates to 
the type o£ features present on the surface of the 
earth. Urban buildings, lakes, rivers, vegetation 
and roads are all examples of land cover types. The 
land use is man’s activities on and in relation to a 
specific piece of land. Land use has been studied 
from many diverse viewpoints so that no one single 
definition is really appropriate in all different 
contexts. It is possible, fear example, to look at 
land use from land capability point of view by eva- 
luating the land in relation to various natural charac- 
teristics namely climate, geology, soil, hydrology, 
topography and biology. 

Since land cover includes physical structures 
built by human being, biotic phenomena such as natural 
vegetation, agricultural crops and animal life and any 
other ty|>e of development, based on the observation of 



land cover as proxy, one hopes to Infer human activities 
and land use* However, there are human activities that 
may not be directly related to the type of land cover, 
such as recreational activities* Other problems include 
multiple use which may occurj^ simultaneously or alternately# 
vertical arrangement of uses and the minimum area size of 
mapping* Therefore, land use and land cover mapping require 
some arbitrary decisions to be made and the resultant maps 
inevitably contain scxne degree of generalized information 
according to their scales and purposes of applications. 

The use of panchromatic, medium - scale aerial 
photographs to map land use has been an accepted practise 
since 1940s. More recently, small scale aerial photographs 
and satellite Images have been utilized for land use/land 
cover mapping of large areas* The development of improved 
sensors to provide high spatial resolution (LANDSAT TM, 

SPOT and future IRS) has made it possible to derive more 
and more information on land cover /land use from remotely 
sensed data* The availability of satellite image data in 
different forms has attracted the flood of applications 
in the land use and land cover mapping field. The advan- 
tage of satellite data are numerous* For the purpose of 
land use m 2 g)ping, the wide and repetitive coverage afforded 



by the satellite platforms are specially important with 
regard to the cost effectiveness of collecting and the 
ease of updating the land use data. Initially, the appli- 
cations tended to concentrate on a manual approach of 
extracting the land use data and producing maps. Witii the 
more general availability of software packages for main- 
frame computer and specialized computer based image analysis 
systems, digital analysis of the satellite data, particu- 
larly those from landsat, has beccxne standard in recent 
years, although the role played by the manual approach 
continues to be important in developing countries. 

1.2 DATA ACQUISITION FOR STUDY AREA 

Ihe study area was chosen as a portion of land 
enclosed between 25°l5'N (lat), 25“35'N (lat), SB^OO'E 
(long.) andSS^^lS'N (long.). It covered bulk of Varanasi 
district and parts of Jaunpur and Ghazipur districts. 

Ihe relevant data for land use mapping of the study 
area were derived from three sources, namely a) the CCT 
containing the brightness digital number (DN) values, 
b) the topographical sheets for concerned area and last 
though not the least c) the 'ground truth'. 
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‘ Ihe computer compatible tape for the area was acquired 
from National Remote Sensing Agency (NRSA), Hyderabad. It 
contained MSS data from the lANDSAT satellite, received under 
cooperation agreements with United States at the Shadnagar 
earth station of I'IRSA, Hyderabad. This CCT covering the 
scene identified by path number - row number 142 - 042 was 
created in four bands inter leaved BIL format on June 8tli, 
1986. 

The toposheets required for the present work were 
made available by the Survey of India. 

The field v;ork for collection of ground truths was 
conducted in the mentioned study area. It vrould not be 
wrong to presume that the land use pattern have undergone 
a change since creation of imagery, yet it was most appro- 
priate to rely upon as the change was not considerable in 
a span of less than two years at least in built up areas, 
land under cultivation, highways and railways. The details 
of field work will be presented at appropriate place. 

J, 3 organization of work 

The course of present work is well described in 
following steps — Step 1 . The number of land use classes 
was decided after a careful study of the toposheets and 
other information available for the area. Once the possible 
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classes tliat could be identified was decided talking Into 
consideration the zresolution of f®s data« 40 trainiiig 
points were selected in each area representing different 
land use classes. The reflectance ON values for thesMti 
training points in all the four bands constituted the 
training data set to be used for classif Icatlcas purpose 
in following steps* 

step 2 m statistical analysis of training <3tota set was 
perfcrn^ and divergences were calculated to show the 
separatibili^ of different classes (land use)* Ttii.a 
divergence test gives a very good picture of the land 
use classes that could be classified tmiquely* It also 
gives the inf cnrmati«xi about the most eccsiapaical band 
cc«ibiiiation for classificatic»i purpose* 

Step 3 •» Eased <m tl^ statistical data cibtainsd in 
step different pixels represented by m values are 
classified as one of the many bmd tme classes decided 
already in step t* A line printer land *> use map is 
prepared once all the pixels are classified* Area 
failing under differ«ait class Is calculated at the 
same time to make obvious the land ues pattern adopted 
in the study area* 



8^m 4 *■» fiaid work gave an oppoctunlty to asses the 
accisracy o£ prepared land us© map* llhe f sets on tlie map 
ie the class undear idiich a particular area was envisao’ed 
to 1?©^ was cross checked with the gromad truth during the 
field work# 
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CHAPTER 2 


PRIHCIPI^ES AND METmDS ZNVOiarED 


2«1 Remote Sensing 

R^note sudsing is the science a£ obtaining infcam* 
tion about an object# area^ cr phcsdcmiadcsi through the 
anaijsis of <3ata acguired by a device that is tsot in c^cm* 
tact with the object# area# or {heiiomeiion under investiga- 
tion* i!he ptiiae objective o£ i^wmote s«dsiiig is to eietraot 
enviromRental and natural resourci^ <iita related to our 
earth* Xnformaticxi about the object cemcemed is conveyed 
to the observer through eleetrooaagnetic energy vhich is 
the infotiKatiesi oarrJUnr and jirovides coiimunication lii^.# 
Thi^ remote sensing data is basloaliy wavelength «* inten- 
sity inCcnnnation idiioh needs to he decoded before message 
can be fully iindtaimt^Mad* This decoding process is analogus 
to the interpretation of remotely sensed imagery idtioh 
ia^pinges heavily on ones- knowledge of properties of IM 
radiati.on* 

the principles involved In em radiatioiis are too 
well knonm to be elaborated here* 0m point that needs 
to be emibasisedi here is that the inteacacticm of these 
radiations with atmospbire aiid earth surface .features 
slKsild be propsrly tsiderstood before hand in order to 
grasp the eesense of remote winsiif,* 
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f'igyre C 2»1 ) ^cshemaLticaXXy iXXustrates the genera* 
XJUsed processes and eXen»nts ixnroX^ed In eXectrcBiagnetic ' 
reflote sensing a£ earth resources* 

Figure (2*2 ) gives overview o£ how remoteXif sensed 
data are turned into useful information* 

H@soXutia:i is a very iJ^xartant x>dr«atieter in resaote 
sensing* Resolution or resolving power is a iseasitre o£ 
the ability of an csptical systeoa to distingui#i between 
signals that are spatially near car spectrally similar 
(Swain and I3avis« 1978)* The ability to laeasure a Mo* 
libysical vacflable using reRK>te s«»sifig requires careful 
consideration csf f omr types oE resolutions * 

1* Spectral* 

2« Spatial* 

3* tesqporal and 
4* radiometric resolin^ion* 

Spectral resolutioei reftce to the dljnension and 
ntmber of specific wavelmagth internals in the electrO' 
leagnetie spectrum to idiich a sensor is sensitive* 

Spatial resolution is a measure of the smallest angular 
or liimar ae^>aratioQ beti^tmi two cbjects that can be 
resolved by sensor* Teiiporal resoltitioKi ^ a; sensor 
system refers to how often a given sensor obtains ima* 




PIgure 2 » 1 Electromaj^nctic remote sensing of earth resources. 
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Figure 2-2 Overview of how remotely sensed data are turned into useful information. Once the 
remotely sensed data are collected, both on-board and ground preprocessing is performed to remove 
geometric and radiometric distortions. This may involve the use of ephemeris and/or ancillary 
data. The data are then ready for visual or digital analysis to extract meaningful information. 
This thematic information is then distributed and hopefully used. (After Landgrebe, 1983.) 



Figure 2- 3 Chronological launch and retirement history of the Landsat scries of satellites 
(1, 2, 3, 4, and 5) from 1972 to 1984. 
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geiry o£ a particular ar«a* RadHoaaetrle riwoltitlon defines 
the sensitiirity of a detector to diff«cene«ai in sigma 1 
str®£igtli as it records the radiant flitr- reflected or emitted 
from the terrain* Hadiaiit flux reccardod hy thematic mapper 
sensors on iandsaTs 4 and 5 x^orded data in S bits (value 
from O to 255) in -six of seven bands# Thus the lather 
system has isptoved radiometric resoluticm* 

Table ( 2#1 ) 0 ires sexae remote sensor systems £unc» 
ticking f rcxa iJanuary 1 # 19^ to Deceoiber 31, 1984 and tM.r 
capability for earth resource mapping*' 

2.1.1 lAKimT f©IinSP!BCTRAL SCAMNER AND THEMATIC MAFPSl 

SENSOR* 


In 1967|, the National Aeronautics mA Space Adteinis* 
tration (NASA) encouraged by os Departtaent of Interliae, 
initiated the Earth Rescxtrcim Technologry Satellite (IRTS) 
prograimie* This programme resulted in the deployment of 
five satellites carrying « variety of remote sensing systems 
designed primarily to acguire earth reuiource information# 

The ^ronological launch subsequent retirement history 
.■«NI some <X! the satellites are ehovn in>ig* C 2*3 
ERTS *• 1 satellite leunshed on duly 23, If 72, vas the first 
experimental system designed to test tbs feasibility of 
collecting earth resouroimt dat# by ummnned satellites* 
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Table 2*2 * Cowiparlscm d laiadsat MSS versm 
Aircraft Areal Coveri^e* 


Platfom 

Scale 


of 

images or 

Landsat MSS 1,2* and 3 

Mia 


1 


185 X 178km »32,930 km^ 
115 X lllral »i2#765 mi^ 


h^tj. 



low-altitude aircraft 

1*15*000 


iooo 

Fo""* ^ 

lotss-altitude aircraft 

1*30*000 


1500 

%/6^ 

Hl^v-altltude aircraft 

1*60*000 


300 

■{Vo 

Coosnerclal Jet aircraf t 

1*90*000 


150 

O 

Govemm^t clviHAn aircraft Is 120*000 


85 

|.^cr 

Gav&crmmnt military aircraft li 250,000 


30 

%1 
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Prloar to launch of ER*rs - 3 on 22, 1975, msh 

rmmmed ERTS program landsat, dietlnguiUilng it from 


the seasat oceanographic satellite launched June 26, 
1978# At this time ERTS «• 1 was ^^troactivel:^ named 
Landsat «» 1 and ERTS •• 3 became lendsat •» 2 at laumsh^ 


handset 3 was launched ffereh 5, 1978, handset 4 on 
July 16, 1982, and handset 5 on March 1, 1984# A 


variety of mechanical failures proi^ted the retlremofit 


of some the landsat satellites# 


2*1*2 SENSOR SySTlMS m BOARD THE lANDSATS* 

S(esisors carried aloft by the landsat satellites 
included the multispectral SK^tnner (MSS), return beam 
vidlcon camera (RBV) and the thematic mapper {TM)» 

MSS, placed on all the five landsat satellites, has 
collected mcore useful earth resource remote sensor data 
of the world than any other sensor* Although the measure#^ 
ment of landscape brightmass is made from a 6241 Sqm 
{?9m X 79m} area, 'sach pixel is sfeformattsd as if the . 
raeasur«HMMit were nailte' from' 4424 Bgpa <S6m x Wm) area 
Cf'ig C 2*4 )}* 

landsat MSB products have provided an umpreoe** 
dmstei ability to observe large geographic area within 




SAMPLING INTERVAL OF MSS 

9.95 9.95 9.95 



Figure 2-4 Relationship between the original 79 X 79 m IFOV of the Landsat MSS and the 
rate at which it was resampled (i.c.» cvcr>’ 9.95 p-s). This resulted in picture elements (pixels) that 
were 56 X 79 m in dimension on tapes purchased from the EROS Data Center at Sioux Falls, 
S.D, 



79 m 
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a single; iinag©* Hiis allows regicsial texrrain anslifais 
to toe periorroed using one data sou3?ce rath®t than a 
multitude of aerial photographs (Table ( 2*2 ) )• 

The theanatic roaiper (TM) seaisor system is a scaniw 
ing optical * mechanical sensor syst^ that reecucds 
reflected and emitted ^energy in the visible, reflective » 
infrared# middle •» infrared# auad thermal infrared regions 
of the electrcsaagnetic spectrum* It collects imiltispectral 
imagery that has higher spatial# spectral and radicmietric 
resolution than the landsat MSS* The plateform for tm is 
stusm in f ig C 2*S}« 

2*2 imoE mocms&xm s^stums* 

Once the remotely sens^ data are in digital format# 
it is possible to analyse them using a digital isege proeeN> 
ssing systfiA idiich can hc^fully extract meaningful infer* 
mation from the imagery* 

Three appeoadhes are often used to p^ovm digital 
Image processing in edncaticMoal amid research ^vircasaants 
iiSmamit 1983 1* The firit is the maiitframe approadh# idsere - 
all analysts utilise a mainframe coRputer (1 32 « hits) in 
batch or interactive mode* They usually vatk with alphanu* 
mmM mmepdM car line printer emtpst aad only occasionally 
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haw the oppcart unity to view and iiit«ttactlvely analyi 
the remote sensor data on a high 'resolntlon hJaelN-.aa 
whit® or coicar mmttox* 

ihe aec<Bid approach involves the use oC a mini* 
ccw^ut^ il6 •• or 32 •» hits) an ine^e procesis«r# and 
high *« resolntioxi oolour loonltor which provide an id 
ijaage jarocessisig environment* 

Finally* theare exist miarocora|mt«ri»toased Il6 -* 
or 32 m hits) i m a g e processing systams* Hiweroixs 
relatively inexpensive microcoE^put^ * based image 
processing systems may now be aoguired £or use as, 
multiple worlestatiofis in a digital image processing 

2*2*1 BRlOimiESS mp* 

' Answer to the Question • do interpreters 
convert the brightness values stcred on a CCT into 
an image that begins to^ approidbmite the mantixmmm 
tone photographs they are used to interpreting 7 is 
the creati<xi o€. Brightness map* also cdomonly rei«rr< 
to as gray «• aisale map# 

%e brightness map is a cosput^ graphic dlspl 
o£ the brightness values iound in digital remote s ew - 
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^ata Ideally^ thera is a mm «» to i« sue 

relationship bi^ween input brightness mlues and the 
resultant inti^sities of output trlghtfiess values on 
display* All Inrightness values from O to 255 wmild be 
displayed as a ccmtlnuum of grays from black to white* 

F<»r pr^paraticsi &i remote sensing brightness map 
two fundamentally differ^t output devicws i line ‘pcin* 
ter hard copy and vidao displays have been used fcr 
present work* Slg {2*i ^ j^ows a porticai of line prin- 
ter hard - copy* 

2*2*2 imss ErmuaSMEHT. 

liaage ehhavzesaent algorithsm are applied to 
reasaotely sensed data to impeove the appearance of an 
iimage for human visual analysis* or occasionally for 
subsequent machine analyois# l^ifferent image eidiance- 
mnt tecttniqiass that were ummS are given be'loif i 

A} Image Hsgnifieaticsm* 

Digital image magfiifleatien (often referred to 
as soomirg') is 'usually |N#rformed to isgsrove the scale 
e£ a display. for visual interpretation purposes* or 
-ocsaasionally* to mal^ the scale of #aotib«r image* 
Cobaan replioatioii repres-^ts the simplest form of 
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digital iwge iia^fni£ieatlc«a» To aaegiilJlf a digital iiaage 
by a factor nt « igaci) piseal in origiiial ixoage is usually 
replaced by an in ac m block o£ pisstlsi. all with the same 
brigbtn<^s wlue as original inpxt pisel* Bm mmeiplm 
of the logic of a four «> tiines nagnificatioti is Musm 


in rig. Ktml )• 

Original image 
0 13 2 

2 4 2 1 

7 8 5 4 

4 8 8 7 


fig *.2*7 

B) Contrast iki^anceiiieiit» 


magnified image 


0 0 11 

0 O 1 1 

2 2 4 4 

2 2 4 4 

7 7 8 8 

7 7 8 8 

4 4 8 9 

4 4 f ' f 


3 3 2 2 

3 3 2 2 

2 2 11 

2 2 11 

5 5 4 4 

5 § 4 4 

8 8 7 7 

8 8 7 7 


i^note seaaors record ref looted and ewltted radimt 
ilmt from earth surfaee materialsNi Tdeally# am material 
eould reflect a trimaendous amount of energy In certalji 
wavelengths » idiile another aaaturial would reflect much 
less energy in the same wavelengths# this would result 
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Aij ec»iti?ast Isetwe^ the tsio typ«>s odE rssaterials idieii recorta«d 
by a rmote causing systems# UnfortmafceJy# different ®at©- 
rials often reflect similar anioyiits amdletit flnx thrcnigiiJ* 
nut the visible and near ** infrared porticms of elastrcfflas^netld 
spectrum# resnltlng in a relatively low contrast image# la 
addition# besides this o3^c«js low - contrast charact^istic 
of biophysical materials# there are cultural factors at worlc# 
For easaample# people in devslopiag countries often use natural 
building materials (eg# wood# soil* etc# ) in constnicticn of 
urban areas# This results in much lower contrast remotely 
sensed imagery -Itxt such areas than for urban areas in develop- 
ed countries# where comsrete# ai^^alt and fertilised green 
veg (Station may be more prevelant# Thus digital image ccmtrast 
enhancement is sev^ly required# Some methods of acccm^li^iJUig 
^ this are otetailed below $ 

8} Slicing and Cray - Scale revise 1# 

Qray scale adjustsumt irules also be specified 
algdhraically# dne grey scale adjustment# in' such comon 
use that it has been named slicing# can be specified as 
follows I 

® * O « 0# I# 2# ..## IC 
w a I « «# lE+l# K+2##### M 
•I d f ■ d « iif^i# # * « t # 
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111# u»«r the threahoMs K and M* G is the input 

gray . level and i* ©utpuft gray level* 

Gray - level ren^raal can also be erpcessed in 
a progfraa algebraically as follow® t 


'adjusted 


** 2S5 *» G 


where it ha© hmm assuraed that both the input aaid output 
gray levels range frcsn O to. 2S5 <8 bit coding}* 


Mhear..S.t:ge.tcai 

Another special form of gray scale adjustroent# 
called stretciiing# is erpreased in a program algebcaioally# 
as follows t 

°aajuirtsaa “ ° * O mOt I, 2, .... K 

• (A ♦ G3 ^ B I O #» m,.**** M 

m O f O » fl|-l|i «•*« 255 
^e paraiaeters A and b are set 
A It 



B • muK * A 

Again;,, this assumes an output gray level range o£ O 2SS 
I 8 «« bit coding)* Ibis linear stretching allows a tligital 






Vi* m. 2*i', 

HKU;HTNK'i* \ > 


CHARLESTON. SC 
THEMATIC M APPLE' 

BAND 4 (MLN MAX STRETCH) 



CHARLESTON. SC 

THEMATIC MAPPER 

BAND 4 (PERCENTAGE STRETCH) 


Figure 2 ♦ 8(a) The histogram of the original band 4 Charleston, S.C., thematic mapper data; 
(b) the histogram after a minimum-maximum contrast stretch has been applied to the data; (c) 
the histogram after a 5% linear contrast stretch. 



iC>udnbay 


iCDuanbaj^ 
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Ijraage o£ restricted gray level, range to tiee ^e full 
dynamic range o£ display device* 

Image analyst® often specify K and M that lie a 
certain parentage of pixels 'into* the tails of the 
histogram* This is called a ■ percentage linear ccwitrast 
stretch* imlilee the earlier one ehich is called mirynasc 
contrast stretch* 

The application of these linear stretches to the 
charleston hand 4 TH data resulted into histograns that 
are shoim in Fig { 2*8^» logic o£ such a linear 
stretxsh is shown diagracmatJcally in Fig < 2*9 )# 
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CHAPTER 3 


TMFORMATIOH EXTRACTION 


3# I Iiati?o<luctJLc«i * 

It is possible to analirsse reraotely s«cised <Sat» 
of tJie ©artb and eactraet useful thematic informatics • 
It is to be noted that data are transfacnaed into 
iirEorraation* cmm of the most oftasi used methods of 
information esKtraetics is imiltis^ctral classifica- 
tiesi* This peioacedbxre assumes that imagery of a speci- 
fic gi^raphic area is collected in oniltiple regions 
of electromagiMttic <^)ectriiia and that the images are 
in good registration (Fig ( j) )• Most of the infer- 
maticsi extraction techniques rely m analysis of the 
spectral reflectance pri^rties of such imagery isid 
employ special algorithms designed to perform various 
types of “spectral analysis"# The procsess of multi- 
spectral classificatiesi may be perfommad using either 
of two methods i supervised or unsupsrvisedt former 
being used for preseat worb* But before we start with 
description of this method^ initial statistics extrac- 
tion of training data set must be dealt with* 

3#.2 Initial Statistics Extraction i 

The mean# standard deviation and variance are 
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tisefiil Bieftstare# a€ eeratral tendency# ifcwreveif# «?fen .fflore 
uee£iil. frcm e ree^te sensing stmdpoiiit &m mmmmm oi 
eovsrlence and cnrreiation anong tlie seyeraX bands# 
flatheiaatlcalXy , 

Standard deyiaticsi in ref Xectance vaXnes in 
band b. 

coyarie^e betwe^ bai»ls k aix3 ckE sas^Xe data 
n 

‘=‘*'lcc “ n lil '®^1X * -*‘k *'*0 ' 

ccxrreXaticsi coefficient betimen the data in 


bands k and c 



I 


Cor, 




®k * ^ 



. they sail pnoride insight intn (M'ta redin^Uuiey or 
guaXity tie^ are the data any ralUie for the reroote 
sensing task at h«aid >• Hi^ coefficient eg cortreXati^a 
tends' to ^iher substantia X redundant spectral infersstion 
in the concertMid 'thanneXs# 

3*3 Supervised CXassificatieii i 

^e foXXouing are ispoctant aspects o£ conducting 
a rigorous and hc^fuXXy useful sup^ised classification " 
of rsBiote sensor i 
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I* to classlficatioKi sdbeiiia imist 

be ^3k»ptedi« 

t* R«pc«se*itat±ire training sites must be selected# 

3* Statistics must be extracted fron the training 
site spectral data* 

4* The statistics are analysed to select the appro* 
priate feature (bands) to be used in the classificaticm 
process# this may Imrolee both eoeputer graphic and/or 
statistical methods of evaluating the degree of between* 
class separability* 

5* Select appropriate classification algorithm* 

6* Classify the imagery into » classes* 

7« Statistically evaluate the classificati'Oii' 
accuracy,*. 

It is instructive to review few of these parameters indlvi*' 
dually to ehow their significance* 

^»3*i * 

Those sites within the image that are reirmsenta* 
tive of the Isaal cover classes of int<mfest may be selected* 
ThesMi sites «3i^ld not be atypical hat ones that r<iq^»resent 
the norm for each class* the image coordinates of these 
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slt«s are then i<lentlflied and to oxtraot statistloa 
frcm the spectral data for each c^ these areas* l^alnlng 
data should he of value if the eiivircmsent froa iitiiGh th^ 
were obtained is relatively hcxaog^eous* 

Each site is usually C€vo^)osed of many piasels* Hie 
general rule is that if trainii^l data are being extraoted 
froia n bands* the ndniiami nuiiier of pixels in a olass Should 
be n*!-!* this condition will allow the inverse of eovariaxce 
aiatrix: for each class to be cabrti.lated, which is in^crtant 
for few classification algorithm* i^ally* > 10 n pinebs 
of training data are collected for each class* 


If four bands are used for inforeatien extraction * 
then each pia^l'in each training ' site is represented by 
a snasuresoent vector# such that 



where ^ brightness value for the li»J)th pixel 

in band h* the brightness value for each pisral in each 
b^ ill: each training class can 'then be analysed statisti* 
csally to- yield a ineMi ^laeasur^aasnt vector# for each 



tab ere "ok rei^eeents the meaxi value of the (3a ta cS^taiueii 
for class, c in band }c» the wm measur^neiit vector c^ 
also be anal]fs@d to yield the covariance vatrix for each 
class c s 


'cJcl 




<^C11 

=Wcl2 

=“%lW 

="<!U 


COW-H’) 

®®*e23 

=^C24 


“'*03* 

C«c33 

°®’o34 

'^*41 

=‘^042 

'^"e43 

^^c44 

is the coveriaiwse 

of class c between bands 


where 

k thrmigh i* J^or brevity the rotaticvi for the covariandii 
matrix for class c (ie \lEl^ will be shcartened to 


^•^*2 gwttlire Selecticc t : 

h judganaent amist be made to determine those bands 
that are most effective in discriminating eaeii class frcM 
all others* this process is conmonif’ ciiled feature seleo* 
tic^f ^ the goal is to dels^ froii the analysis those hmdm 
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that prc3«ride only iredtja<3aiit spectral ii5f€3se»tic«i* In tills 
way tlie disaensicmality (le# the nuraiber csf bands to be paro* 
cessed) Iti the data set m&f be r@d-ij«Ns3« 

Statistical i^thods of feature selecticKi^ 'mem used to 
qumititatlvely select the subset o£ bands (or features) that 
provides the greatest degree of separability betwean aatiy two 
classes# c and d* Diwergenee was one caP the first aieaswre ’ 
of statistical separability used in the machine processing 
of rsdsote sensor data# and it is still widely used as a 
a^thod of feature selection* It addresses the basic problaii 
of deciding what is the best <|iwband siibset for use in the 
supervised classificaticm process* Ihe ninifber d£ 

tions# c# of n b«sids tahen q at a tine is 

# n % n I 

dm) m 

^ ql (ii«»q)f 

fhe degree of divi^igence or "separabilitf'* between 
two classwi m and d# Qiverg^ # is computed according to 
the fonaula 

niv«g^ - |*T 

' T " ■ ■■ '' 

+ 0.5 * 

%itiere {*i ’ is the trace of a inatrix Cie* the aim 

w 
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o£ cSiaigoiial. elements and are tdie covariance matrix 
fco: two classes# c and d# tinder investigation# and ^ and 
t5^ are mean vectors for classes c and d* 

If there are more than two classes# the average 
divergence, I3iv@cg _^ , . , Is coisipnted' over all p<»islbie pairs 
of classes e and d# while holding the stftwet oi hands# g 
c^^stant# as follows a 

I m-*l m. 

Dlverg*^ » * £ % Etwerf 

avg c d*c+l 

Komar and Silya <1977) enggest that It Is possible 
to take the divergence logic one step further and coi^te 
transfoacmed divergence# expressed as 

T C^Edverg^^) 

Ddverg^jj^ m 2000 £ 1 • eaqpl * - | ] 

Ibis statistic gives an esKponentlallir decreasing weight to 
Increasing dlstaiKJe 'between the class# It also scales the 
divergence value to He between O and 2000» 

h transfcarmed dlveegenc^ value of 2000 suggests 
excellent between -* class separatloiit jitoove 1900 provldahi 
good separation# while below 1700 la poor# Kor divergence 
values less than ISOO we get spectrally similar classes# 
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3»3*3 Claissificatioa a.l^ogittiai 5 

VaarAous supervised classAficatAou' Hsetbodls raay to© 
used to assAgn ©u isalaaaHU pixel to o«e <£ © lawtoer o£ 
classes* The choice of a particular classifier or deel» 
slou rule depeuds m the native of the Input data and 
the desired output (i^iedraan# 1980)* Parametric classl-* 
f icatim algorithms aasiviie that the otoserved measuremefit 
vectars 3^ c^^lnsd for each class in each spectral toand 
during the training phase of supervised classification 
are Gaussian In nature tl#« they are nmaallY dlstrltouted)* 
Hcnparametrlc classification algorithms ma3oi no sirh assunp* 
tlon* For present vorto maxlmm lileelyhc^ decision rule 
'whl^dbt Is a special case of Bayes* declsloii rule* Is used* 

ihls decision rule assigns eadh pixel having pattern 
mmaniBcmimtm or .imtmm X to> the class 0 ihose units are 
OKset larobable esc Hkely to have given rise to feature 
vectea: X* 

It can toe provad mathematically that this rule 
minimise the **^Ktotal %ror oi Classiflcatloi}*** If there 
are *G* groups cn: classes then according to Bayes* Ihile* 
assign the 4Tieasurenent X to group t t 

)► FCGjj/SH) for : all 1 J* 



Hcsw according to “Bayos TtitPcrcBa of Conditioned ProbaMUty** 


PCX/S.) P(G.) PCX/S.) P(G.> 

. gssr-.^ an -«V<5j!Sf5p - ^ p{^, 

wti'ere ** l^oibaMlity of cls.»® '*1* 

*l1ius the final form of Bayes Rule Is 
R{3{yt?|^) RCg^J > l?{o^5 

for all J 1 


To apply Bayes * Rule ym have to Jsnow the yalne of aeid 

for ©acti group* It can toe i^am that it Is easy to 
find the vail*©'' of *Mjt so far 

How If Me assmae that if within each group the variatolea 
that mahe np the measuremaat venter X have a multivariate Horm 
dhatritoiiticsi# then ihe form of ccwKlltioiial sardbetoility 


is given toy the density function a 

1. 


P(X/0|,) m 


* , , -'I ■ t 


|vj 

Wh'ere m mean of class of groip 1 

m varlanoe covariance iwitriac of class *1 
p(Gjl w A priori ppofeability of class *i* 

■ . ...I 

(Xm^LFm TrastSpoee of 

'St ■ 
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•aerofore estimating Vlx/0^) oa»s dc«i to estimating 

two pairaiMsteor®# tiie groij^ siean vectcac and the giroias> 

covariasice matrix. t£ there are n fasasurements* in X 
or channels) then there are *»* s^ans in the i^aan 

wsctcr «rK3 eieraents is the ccatariance matrix making 

2 2 

a total o£ quantities to be ^itiwated for each group 

tout ^is is coneldterably fewer than the »u«i>er of requiri^' 
for a direct application of Bayes rule usio^ estimates of 

y 

relative freqtjsnoies* 

Using the normal foami of 

have t 

Assign X to Gj if 

In |Ti| 4(X«Mj,) iXm 

i 

j In 1^ 1‘i' iXm 

for all J i 

Although matrix operatlcm are involved in working mt the 
atoove exxxreesionsi^ the final result is a nuslxst fiar eadh 
group Which depends m X« Ibis quantity is c^tm referred 
to as .ajscriiainiest. ....^Oige. and the expression % 

4^iX) m S»|?J j|^) (Xfc i|^)^ 


p(X/t5jl^) in Bayes rule we 


t 

M|) ^2in P«0|^)| / 

■U j (3fi» »2ln P(Oj) 
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is called a ‘• Mscriaiiiiant giaicfcicm ” * 

..£ig?a.l rcle for Bayesian mettoCKi classification^ 
J^ssign the laeasur^nent X to the groij^ with the siaalleBlp wmlmt 
of PiQ^)p we can infer that as gets SMail®atf» 

the evideiKje of raeratoars of gronE> *i* inoareasea %daidh is the 
opposite to the hehavioor of PtOj^^/JC)* Syen though does 

not b^av© lilce Pto^/JC)* It is is^Kartant to srealise that the 
rnle » •>• 

Assign to groiqp ‘i* if t 

d^(x) - in p{o^) ^ d^{x) ^ In piop for all J i 

is the Baye's rule if X is .BRiltirariate Homally distrihuted 


in each of the classes or g£«n^# 

ihis Bayes* decision rule is identical to the isnsciratum 
lilcellhood decision rule msm&t that it does not assuma that 
«a«ii of the class is equally x^cfcahle* A yrlcari proibahilltles 
hare been need sucsejesfuily as a way of incocporatii^ the effect 
of reledf , and other learrain dharacteristics in laprorlng 
classification, cKxnirecy* 

this type of error analysis eos$«res the scseuracy of ' 1 

the remote S€S)sinf «• deriired classifiectioii nap plseal by piKei 



39 


with tli$ asstitned tr«a laud aaap* h ninS^er of strategies 
have h@m develc^d# 

first# it is pt^sible to ooiiduct a site • specific 
erarcar evaluaticsi based only m the training pixels need 
to train the classifier in a supecvised classlficaticsi* 

This singly means that those pixel locations i»J tised to 
train the classifiers are carefully evaluated on both the 
remote sensing » derived classif Ication map and the true 
map. If training saiaples are distributed randomly through-* 
out the study area# this evaluation may be considered 
representative of the study tftifiMrt«aitely# the 

locations <si training sites aie usually .nonrandon *•- they 
are biased by the analyst’s a priori knowledge of where 
certain iaiid cover , types escist, in the scrnie# Because of 
this bias# tlie classification ac;curii.cies for ptasels fowajd 
within the training sites are g^^ally higher tl^an fc» the 
remainder of the «fap because these are data locations that 
were used to ** train** the classif<3r» 

Conversely# if other test locations in the study 
area are identified and correctly Istieled prior to tlie 
classification and if these are not used in training a£ 
the classif ication algorithm# they can be used to evaluate 
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th€! accuracy classif icaticn mmp* Hiia procoSurc gniiew 

rally yields a lac^ credible dasalflcatlou accuracy assess- 
sient* Ucimwer,. additi<^3, **groia 0 ^ trutih'* is req*^®<3 for 
these test site areas# coupled with the prdhl®® o£ identi- 
fying how many pjbeela are nec^sary in each t<^t .site class • 

Ihe ideal nunh^ c£ points to be t^ted in the land - 
use classificationi map was determined from the formulas for 
the binoraiai probability theory* Hie fcaraiila for the niaii^er 
of points selected was s 

H « I i 

Where p is the expected percent euccuracy# q the difference 
between 100 and p# B the allowable «ra!rrimr# and H the nusSaer 
of .points to be sallied (Fitspaiuricle - lt@o)« 

Once the criterlou for 0b5«®tiv«ly identifying the 
location 0 ^ specific plrels to Iwe compared is 'detersdned#. \ 
,it is necessary to Idiwtify the class assigned to. each data 
hre talmlated and aeepcrted in a contingency table# niiere 
overall classification acrjurascy and raise lassifieaticai 
betwe^ categories are identified* 

are .two types. of errors t 

I* ■ "li pdbfiol may be assigned tyo a class to which it 
«^es moft belong C«ffi'error of oewaission)# 



a» A piaml is assigfned its af^ropfiato 

class Caa mxrm: of omiaaipi)* 

Th« a£f«»diagofial. entariaa prwida info»®ation on 
oinrora of wsission and ccwBiisaion#, Errors e£ ccawaisslon 
ara found in ttia lower laft naif of tJie laatiriat;# tiftiila 
errcrs of onidssioii ara f tmnd in tii® ui^r, ri#»t half 
of the natriac* Sferors of oailssion for eacti elass ar® 
coRiputed hy sunffirfjEig . the norafoer of jji^cels assigFfieai to 
Inoorreot categories along each ra? and dl’ifdding this 
nuntser by th© total nvriber of trtie pfccels in the category# 
Errors of casosissioai are ccxitpited by siwxitig the m®i»er 
of piscels assigned to incorreot categories along each 
eohma md disiding this nn^b^ by the total nnolaer of 
piaeebi assigned to the ool»n category * 
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CtmPTER 4 

CIASSaCgXCAliQK SCHEMES 


Qae crucial factor in detetiaining the smsc&ss of land 
use and land cerrer mapping lies in the cSioice of an appro** 
priate classificaticn s<iieHte designed fenr an inteniiid purpose* 
A good classification sdiienne ^ould Ibe easif to use with no 
an^aiguity in defining each land use and land cower categocf* 

It imiet also l&e able to generate the degree of details 
required* Certain classification sch^aes hawe been developed 
that can radily incorporate land use and/or land cover data 
obtained by interpreting remotely sensed data* naiaely US 
Oeological survey land use/land cover classification systeasi, 
the Michigan classification systee# the Ccifardine Wetland 
classification system* etc* 

Major points of difference between various classifica- 
tion scheuMas are their Cl) ei^hasis# and (2) ability to 
inooriKnrate iii£ormati.c^ Obtained using remote seisii^* For 
casaw^le# USOS CAwfersem et al# # 197© ) # is "resource** oriented 
in contrast with various "people or activity" oriented systew 
such as Standard hand Use Coding (SIW) wmial Cihftoan Renewal 
Administratioii* 19^ >* the vms rationale is that "although 
there is an c^ious for an tn^beyn ** ccriented land «• use 
classificatioii system* there is also a need for a resouren *» 



orieiit«€ clA«el.fication syst@fa primainr ■ ©raptoasi# 

wouia be the leeBaaining fS percent of tJnited States land 
area '*• *Ibe syst®n aadresses this need witdi eight of the 
nine level t categcaries treating land area that is not 
in urban or b«ilt«*«p categories (Table 4*1)* The system 
is designed to be (briireii primarily by the interpretatic»* 
o£ remote senscor data obtained at yariotas scales land 
resolutions (Table 4*2) and no>t data collected in situ* 

It was initially develcped to include land use data that 
was visually interpretei# although it has been widely 
used for digital aailtispectral classificaticm as well* 
lend use/land cover data classif ied accordlfig to this 
system is compatible with land use coding used in the 
us Geological Siarvey*s Geographic Infcnimatic^ Retrieval 
and Analysis System (GIRAS) (itLteiiell et al*, IfT?)* 

The USGS l^d use and land cever classification 
system was dssigned according to the following criteriat 
( 1 ) the miniMSR level of interpretation accuracy using 
resKn^e mmmxt data should be at least 85 percent i (2) 
the accuracy interpretation for the several categories 
should be about ©luali (5> repeatable or ^petitive 
results SticmM be obtainable from one interpreter to 
another and from one time of s<misiiig to iMotberi (4) the 
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Te®i3a 4«1 1 U.S* Geological Swcv&^f Land XTse/Lond Co^E* 
Classification System for tise wida ifeaaote 
Seaasoc Dstta* 


lievel I 


jewel II 


i Iflidbaii or lmilt-*tp land 


2 Agricultural land 


3 Hangeland 

4 Pcr^t laiKJ 

5 Water 

s Wetlaisi 
?' Bacren land 

8 Tyi^lre 

9 Fereamial anoir and ice 


11 Kesidcntiai 

12 CcHsaercial and senrices 

13 Industrial 

14 Transportation*, ecwaaaiiicatlon 
and services 

15 Industrial and coscBnercial 
cospleoees 

16 J^tred urban or built-up lanl 

17 Other urban or built-u^ land 

21 Cropland and jasture 

22 Orchards* groves* vineyards* 
nurseries* and omamaital 
bofrticuXttnral areas 

23 Confiiwsd feeding operatlcais 

24 Other agricultural land 

31 Herbaceous rangeland 

32 Shrub and Ircush rang@la* 3 d 

33 Hired rangeland 

41 Qi^ichicms forest land 

42 Evergre^ for^t land 

43 Htsced forest land 

51 Streams and canals 

52 lakes 

5 3 Reservoirs 

54 Bays and estuaries 

61 Fcrtested wetland 

62 !slc«ifc 0 pested wetland 

71 Dry salt flats 

72 Beadhes 

73 Sandy ajreas other than beaches 

74 Bare expcKwsd rocks 

75 Strip mines* gparries* and . i 

gravel pits i 

76 TTansltlcmal ar<M 3 i 

77 Hired barren land 

81 Shr«b and brush tundra | 

32 H«Ebaceous ttmdra 

83 Bare ground 

84 His^ ttaidra 

81 Perennial snowflelds 

82 Glaciers 


Scnirce i Anderson *st 
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Table 4# 2 


•Ebe four levels of tbe u.s. Geological 
E*ana Gse/Iena Cover Classifloatlon System am 
the Typ^ of Remotely Seti^Gt Data Typically 
ifeed to Provide the a^cnaation. 


Classification level 


Typical data characteristics 


lehdsat Cfcormerly EiTS) type of dat^ 
Hlgh-ialtltiaie data aeqiiired at 40 ,oq, 


ft (12«400 a> m: abovei results l® 


Imagery ihat Is less than 1 *00^000 


scale* 


fMliaa»altlt%tde daita acquired betve^ 
10,000 and 40,000 ft (3100 and 12 ^ 4 o,j 
mh results in imagery that is betw 


Is 20,000 and ls®o>ooo aKsalei 


l<a»*altitiide data aeqtiired below lo , 
ft (3100 ro)f results in imagery 
larger than 1*20,000 scale. 


Sources et al# (1976^ and Jensen et al. (1963 }« (q) 

^onearican Society of IhotograaBaetry aaad R^mate Sensing 
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elassifAcaticsa system should lae ai^.liosl>le over ejcteosive 
areasf (S) the eategorJ^atAoii should permit laud use to he 
inferred frtjm the laand cover types# (6) the classification 
system should be suitable for use with reaiote sensor data 
obtained at different times of year# t*J) categories should 
be divisible into more detailed siiflbGategories that can be 
obtain^ from large scale irngmiy or grcaind surveys# CS) 
aggregaticn of categories must be possible# (9) ccm^ariscst 
with future land use and land covecr data ali<n 2 ld be possible# 
and (I0> multiple uses o£ land should be recognised 
possible* 

ihe system is designed to use four "levels" of 
information « two g£ which are detailed in ’Ibble 4«1« 

The I7SGS classificatlisi syst^ also provides fcr the 
inclusion of more detailed land use/land cover categories 
in levels ttt md 1V« fhese two levels can be utilised 
to prcnride ii^omtion at a resolution appropriate for 
regional (multioNnintry)« country or local planning* and 
management activities* ' It $m intended that’ levels ill 
aiB^ W be dasjj^ied by the local us^s of the CBGS system* 
betting in^ mind that the mtegmim in each level must 
aggregate IntO' the ihtegwies in. the nest higher level* 

Fig* (4*1) Chows a sae^le aggregation of classifications 
for levels fix* tX and I* - 
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C<iwair<2in tt alv ^ 19*19} dev«]U>ised a wetland «slasi»i*«- 
fication system that Irascarperates inf ormatiian ■ extraeted 
freei. remote sensor data and in sitn measuremejits:*' ' It 
descrllaes ^ecological tsam^ arrangres 'th^i in a sys tem 
useful to rescaneoe managers, and provides imif<a*aity 
of c<x«2epts and terms CFig (4*2 >)* s«et lands are classi- 
fied based on plant characteristics.# soil and frequency 
of flooding* 

Finally# it Should be noted that there is a 
relationship h&tmmm the level of .detail in a classi- 
fication scheme aiKi spatial resolntioi} of remote s^sor 
systems used to provide infoimation* For wcample# Welch 
(1982) numnarlsed this relatiooi^p for mapping 'fsi mfb&n/ 

Mnd use and land cover in tmited States (^ig 4*1))* 
A sIMlar relationship exists nhen ntaiping v^ititioa (i^'tSAi, 
19i3)«i ' For fchnple# the sensor syptems 'md spatial resolu- 
tions useful for disoriiainating vegdtqition frcsi a glisbal 
to an iU' situ perspsetive mm mwmmri»mA in Figure C4»4 )» 
ibiS' "S'Ug^lswis that the level of details in desired c lassie 
fioation systest dictates tlie spatial resolution of ramote 
mmam data that hould Tm taped* dpeetral resolution is' 
also an important consideration* iiisreiiNHr# it is laot as 



WETLANDS AND DEEP- WATER HABITATS 


SyiUm 


Subnyvttm 


ClaM 


r---MarijrMi 


Subtidal 


—Rock Bottom 

— UnconaobdaUd Bottom 

— Aquatic Bad 
*-R*af 


Intartidal 


■Aquatic Bad 
RmI 

Rocky Short 
•UDConaolidatad Shora 


-Subtidal- 


o EstujLrine • 


- Intar tidal- 


-Tidal- 


-Lowtr Partnnial - 


- Riverine - 


-Upper Partnnial - 


-Intaxnuttani - 


“Rock Bottom 
— Uncon toUda tad Bottom 
-Aquatic Bad 
“RmI 

- Aquatic Bad 
“Ratf 
-Strtambad 
“Rocky Short 
-Unconaoiidatad Short 

- En>arf«Dt Wttland 

- Scrub-Shrub Wttland 
“ Foraatad Wetland 

“ Rock Bottom 

- Unconaoiidatad Bottom 

- Aquatic Bad 
“Rocky Shora 

- Unconaoiidatad Short 

- Emargent Wttland 

“Rock Bottom 

- Unconaoiidatad Bottom 
“Aquatic Bad 
“Rocky Short 
-Unconaoiidatad Short 

- Emargent Wttland 

-Rock Bottom 
-Unconaoiidatad Bottom 
-Aquatic Bad 
“Rocky Shore 
-Unconaoiidatad Shore 

“Straambad 


-Lacuatrine- 


- Limnetic * 


-Littoral - 


£ Rock Bottom 

Unconaoiidatad Bottom 
Aquatic Bad 


f—Rock Bottom 

Unconaoiidatad Bottom 
Aquatic Bad 
}“ Rocky Shore 

Unconaoiidatad Shore 
Emergent Wetland 


L— Palustriae 


—Rock Bottom 
— Unconaolidated Bottom 
—Aquatic Bed 
— UnconadUdatad Shore 
‘ — Moaa“Lichen Wetland 
— Emergent Wetland 
-Scrub-Shrub Wetland 
^ Forested Wetland 


Fig4 #231 The Cowardin ct al. (1979) classification hierarchy of wetlands and deep-water habitats, 
showing systems, subsystems, and classes. The palustrine system does not include dcci>watcr habitats. 
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ITLANI) PALrSTHINK 





PALirSTRlNF nn.AM' 




s TEMPORARILY FLOODED 
b SEASON AULY FLOODED 
c SEMEPERMANENTLY FLOODED 
d INTERMITTENTLY EXPOSED 
e PERMANENTLY FLOODED 


Fig4 • 2t) Distinguishing features and examples of habitats in the riverine system. (Cowardin 
Cl al., 1979.) 
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LAND USE CLAJSS aevel) 


Figure 4 « 4 Spatial resolution (IFOV) re- 
quirements as a function of the mapping re- 
quirements for levels I to IV land-use classes 
in the United States (based on Anderson et 
al., 1976). Levels 1, 11, III, and IV information 
are normally derived from satellite, high-, 
medium-, and low-altitude image data, respec- 
tively. Note the dramatic increase in resolution 
required to map level II classes. (From Welch, 
1982; Jensen et al., 1983.) 0 American Society 
of Photogrammeiry. Used with permission. 
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LE\TEL 1 : Global 

AVURR 

rrmtiutiofi 1.1 km 


LEVEL n : Continental 

AMiRR 

Landsat Multispectral Scanner 
resolution: 1.1 km -- 80 m 


LEVEL m : Biome 

Land5 Multispectral Scanner 
'Hiemt-.c Mapper 
S>Tithetic Aperture Radars 
resi^ution : 80 m — 30 m 


LEVEL IV: Region 

Hiematic Mapper 
High Altitude Aircraft 
Large Format Camera 
SPOT 

resolution : 30 m -- 3 m -h 




Boundary Waters Canoe Area 
/ \ 


LEVEL V: Plot 

High and Low Altitude Aircraft 
resolution : 5 m ^ 1 m-t- 



LEVEL VI : In Situ Sample Site 

Surface Measurements 
and Observ'ations 



Upland Forest Wetland Bum 


Figure 4 «3 Relationship between the level of detail required and the spatial resolution of various 
remote sensing systems for vegetation inventories. (From NASA, 1983; Botkin et ah, 1984.) 
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critical a paraiaetee m spatial rsselwtioR siwse misst 
0f tliwe sensor sfmtem record ^eisgry In approrlraately 
the saiae visible and near infrared portieais oii th« 
electramgnetio speetrnra* 
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CHAPTER 5 

DEVELOPMSHT OF mOCRAM 

S<^%ware dwelc^xment was laost essential laart «f 
present thesis work* It is neeciless to reiterate that 
lifiis work was Imssd on digitizeci data obtained in form 
of C€3®ttfaater ecssK^jatihle tapes (GCT’s) and msiiaal appc^sadte 
f csr data iutocpretation throtigh analysis o£ these data, 
which was enonaotas would ha^e been an unnerving task 
and all tiie b^efits of t^ing satellite platjg^fors® for 
data collection %#ould have been eroded# In idiis cootojct 
the swiftness and acciaracy associated to digital computer 
in land use/land cower mapping ou^t to be appreciated* 

5#1 Software for statistical analysis# 

Statistical analysis of tnrainix^ data was perf^^smed 
by program staisin* it «wed dif f <^ent suiMfoutines to get 
the Maximum, minimum, mean^ standard deviation, cowar iance 
and coefficient ■©€ c«nrrelati<xi between brl^tness value 
= digital nusfoer associated to different pixels in a parti- 
cular class and particular spectral baoid* The functions 
perf cawed by different suhrcnitiiies us^ed in this program 
are listed below* 

STATOl <- parametric inputs to ■'this sdbroutiiic are ni»i>er 
of classes {ICIASS) and nuatser of pixels used for training 


the .classifier to ItlentlfY aasy particular class CitaifA)# 

This siihrootlfie reads the brightness values of training 
points in array fD* it perfcanas the statistical analysis 
and rettims the rainiBura.#, xiaxiittsa, aean,* variance •• c©». 
variance, inverse variance - covariance .matrix and eoef.f- 
icient of correlation matrix*. At the saaie time it returns 
the deterniinant o£ variance • covariance matrix to be used 
for classification stabrontine CSTAitJ3)*. Ibis siibraatine 
also laalces use of two library subroutines (FoiAAf and 
F03AAF), 

STAT02 •• this subroutine perform® two multiplication opera- 
ticais given below 

mo * CRD.l'^ A. |RD|: 

v®2tc3r |RDI and matrix A are input to this stfcroutine 
to mroduee mo as output# 

STAT03 - this subroutine uses the training data to Checlc 
the ef ficlfjncy of classif^lcation algorithm (maxiRaim Mleely- 
hood in .this case) and presents the data in a tabular fcscw 
of ec«?fusi€«i matrix* IwpOiftant inputs as mean# cwariaiK;# • 
matrix and training data itself .are transferred to subroutine 
through labelled cceamjn statement. It classifies the training 
points as belonging to one ssi the classes fremi thefy 

are extracted* 



STAl^ « this siihrouttine jarepares histcigirss®® caf buightniess 
value in a band versus frequency for each class* Siwilta* 
neeusly this subroutine plots the histograms co line printer, 
each *** representing a single data point* 

STAIX© -» giims coincident spectral plots in di££er«it baiwSte 
for different classes* 

STAT06 -* subroutine perfcanas Divergence test for ctotaining 
spectral separability of classes* it tabulates the diver** 
getKjes saod transformed divergences along with their average 
for different b^ad ccnbinatich* This table serves well for 
feature selection* 

Ditfsrent possible band conhdJiations mui class coebina- 
ticns are g«ierated thrcmih subroutine ST&TO?* 

lilnVages of these sihroutines with either sutwroutines 
or main program is shown In Fig C5«lK 

5*2 Software fear classification algorithm s 

This pacl^ge is used for classification of piste is to 
one of the pce-determined classes* The program BCIASS has 
the optltm to use either roaxteum likelyhood algoritha or 
Bayes* algorltlua,, based on a priori probability, if aprlori 
probability for all classes are same the algorithm used is 
the fcarmer one# else the latter caa® is used* Bata input to 
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main progiraa is a vector of reflectance values £roai a part 
of scan Hue in all four spectral It l-s retreiveH 

fro® tlie CCT with help a si:£broutiiie called ailSoi* Ihe 
functions perfor«ed by other subroutines are as follows — 

BCl^l • this subroutine p«rforBis the statistical analysis 
of the training data set reccacded as a ittstriaE After 

the statistical operations the mean vector TDHIAK, ccavariance : 
matrix TDCVAR and inverse covariance matrix: aire available,. 
Since covariance imtrix is not required in tbe iJoforiKiticii 
extracticau process and to econotnize the memory space the 
Invorse covariance matrix is overwritten over tlJ.e existing 
covariatxie matrix TDCyAR* This subroutine also calculates 
the maxlimm eoid minimum of ref lectsmtos values In a partJ^lar 
spectral baixS* All the results of the analysis are transferred 
to the ready subroutines throu^i camixxi stateraent* 

BC1S02 *• this subroutine is the mc^t i ia pcr t ant <sie as it 
contains the classification algorithm* The input to this 
subroutine is as a data matrix R that contaitos the 
-reflectance values of the part &£ 'line that is to 

be classified* The 'probability of a partieulaitr pixel to 
fall within the range of a class(l^O) is calciilated and 
'^ased on this prch^llity the class to which this pixel 
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r^rescfits is dlecided# & cwerstruck -» dbaracterUne •• 
printer land use laap Is generated in tijis subrcmtitoe 
only* snbrcHitlne BGlsos is used in this subroutine 
for raatrix operations* 

BCLS04 - this subroutine transfers the required part 
four coosequentive records ie four bands to subroutine 
B€LS02« The length of r« 2 ord or the nuidser ctf pixels 
to be analysed is giren by variable IlENTH* fhe data 
on CCT nAiich is in binary raodto is changed to deeiraal 
mymimm fear use in analysis*' 

Subroutine BCI 1 S 05 plots the hlstograins to check 
the validity csf the assumption pertaining to Gaussian 
distribution of errors and svSsrxmtSme Wl^36 gives eo* 
spectral plots to i^ow the enpanse c£ ref lectano© values 
in each band for each class* 

Subroutine llitoge is i^own In Fig (S*2)* 

There are other preagrams too that tstere used in 
the present work# iBiey are described below * 

CHPhOT -• this pstegmm has two opticats available*, the 
deeisio® is iBad#' bas^' on th® variable EEICIIII?* If 
DBCnro is •DBH* this program pcarforms dwisity slicing 
based on the mimmtma required asai their tapper' lew^r 
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lindts. If DECIND is *HRD* the sreflectatsce mlye irange 
is divided into sixteen classes and their up^r and IcMer 
Undts are decided hf the ^aro^am* It was used daring 
selectic*! of training data points* "The <xitp«t of this 
program is generated into output files ’s^ere 

n varies from l through 5 depei^ing «an the length of 
record toeing analysed, 

TAPEI'ffl - this program reln?eives a record from CCT 
and changes it to decimal lacxte from the mode of recorv 
ding ie Binary, It tatoulates the reflectai*^ values 
in classes of class interval of 10* Hie output is in 
file FtEEvmT. 

TffTRNS •. this prograaa changes the latitudes and longi- 
tudes of a point on ground into Hue nusdier mis3. pixel 
nimtmr to locate it «ti the scene# It uses csrthcmarpiiic 
projection# 
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CH&PTER 6 

RESUIffS AKD DISCUSSlOKf 

In this chapter the results chtaineS ireia the 
statistical analysis of training data* the classlfi-. 
cation t3£ pixels and the field siarvey a®® presented* 
Discussions ihererer required are also dealt with. 

All the presentable results in foE» of photographs 
and co^uter printouts are iixsluded in this chapter. 

A study of 1*90*000 tbp^heets of the study area 
as discussed in cfeaptea: caie clearly Showed the distinc- 
tive classes as water hcdles* agricultural land* builtup 
area and barren lai^. Few oE ^ese level I classes were 
further pcmsible to be separated into level ll^ For 
eaoifflple* water was reclaaslfied into two level 11 c.la8ses 
based upon the water depth* It was plainied initially to 
separate the Orchards f rooi croplands aie^ pastures tov 
attain level IX classific^ticn of Agricultinral land. 

But the spectral inseparability cf these two level XI 
classes because oE coarse zesolution- of reiaate senscac 
data foiled the atteept. Um&e the classificatiaa' achi- 
eved was as ati'Cnm below t 

level 1 level XX 

1* tlidasn or built-up aar^ , — 

<lncl. roads and railways) 


l.€fyel I 


lievel II 


2 m Agriciiltiital I«ii€ 
5 • Matear 


7* Barren land 


51* Streams of langco: <S^th 
52« lelces and str^eami wi.tii lean 
depth water 

73* Sandy areas lay side e€ the 


rii^ 

77« rni^ed bai^eii not issed 

for agrlcultwe regularly* 

6*1 Statistical Analysis t 

llie statistical analysis results are presented £rc»a 
page to page* fhougti th0 tatbles are mM ea^lai^tory yet 
few points need ^erplailnaiicii liiich will be diseased* 

The correlation coefficient calcmlation hetwean diff** 
erent spectral hands r^aesesiting a class giwes very striking 
results* It nekss note of reteftdancy cf i^peeti^l hands in 
providing further information* Wm eraiiple* in class 1 
(Mater, SICbsgs)) the coepfieteits ^ ccirrslatitni hetween i 
hands seven and fire (0*66) and hands five and four (0*61) 
are helow o,'^ hence the spectral . resiNsnses ,©f both the 
hands in these hand sets do not follow the same pattern* 

It «isures that the different bands produce informations 
that are coH^jl^aentary to- emdi.* It is true fear ^every 
hand for all the classes* 



1h« hlstogira® s^JLysis searves well to clsecl, tSie-, ; 
assui^jtlcti that distrihutlc^ of the clestid poJjatsi^ 
fcHcmlng the categoary tralRlng data Is Gaussian (normally 
distributed)* Hiis asstauptim <aE ncacmalily is g^serally 
reascKiable for ccmon speetral response distribotioii* 

'live histograms for all the bands aare Gaussian in nature* 
though the spread aiid location characterised by earianee 
and mean walu^ of brightness JM are different for each 
histogram* 

Codncidcht speetral plot is another itay to presumt 
the trainiiig data*^ *llie ^read M brightness ealues fi«r ' 
training pcKints in each class are present for s^pectral 

band* Ihe cjyerlap of these spectral plots for different 
classes may hint that the s^^jarability c£ the classes 
whose cospectral plots cjcerlap Is poor in that particnlar 
bai^* For eacai^le* in band 4* the s<^?«cabillty of classes 
Imiltuip area (regpr«s®srted 3) and cultivated area 
(represented a* 4 ) and balltai> area md wat«r of less 
depth (reppes®**^*^ ** 5th class) are questionable* ®imi* 
larly specstral plot for sandy area (repcesented as class 
2 ) and cmltivated area overlap in all the foiir spectral 
plots* 100?S overlap being In band 7^ Burt this plot alone ; 
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can not be the conf lmatory stateiaent about the spectral 
separability. We will haire to ^ for diw^?genee test to 
ensure the separability of classes. 

ihe cJivergence and transformed div'ergresce test is 
conducted fcr two imrposes, namely Cl) separability of 
dif£er®at classes and (2) most cptiwua b^d combination 
for classificatiem algocitlm^ ie feature selection. Hie 
test ccaidueted for the first purpose helped in a way that 
It told in advance that further classificatlcn agrieul- 
tural land in level II classes of cultivated land Crepc®* 
sen ted as class 4) and orchards Crepaotseaattei as class 6) 
is not possible as their diveogeiwse is well be:low the 
rec^mended limit of i&OO* 

Hie featmte seleeticm is the 'Selecti<ni ^of spectral 
bands (features) be used in classification td give best 
results in least of coB®>utati««ml' effort* Hie differ^^t 
band cosbinatioQs (equal to 15 ) are genwated In program 
STATAJt?. For each band cofiabination the divergences vwre 
calculated for all possible class pairs* Hie average 
divergence for each band ccabination was calculated at 
the saH» tiiae. Ihis average divergence gives a measure 
of utility of that particular bs®d cosbination as the 
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feafcujre® in classi£lc^ti<Mi process • the 

of all tJie four bands give the maxtimsB of average txans* 
fcastaed divergetKse bei^;# can Im used witb greatest relia-*^ 
billtjf In class ssparatiGHnu At tbls point it sfeoaM be 
en^phasised that higher niastier of bands used in elasslfi^ 
cation do not always ensure better results in tenas of 
separation of classes* Fcjt example# band 7 aloiie gives 
better separability idten conpared to bands i and 5 ca? 

4 and 6 axe being used in diverg^ice test* 

Etxpressions used for calculations of div^Dgence 
and transfonaed divergence are different and so are the 
acceptable liisits to classify _ two classes as -distinct 
classes* But it seffion the limit 1SO0 for divergenee 
tends to give very -conservative Itroit* Umucm 'transformed 
divergence is vmty appropriate to be used because the liioits 
are in accordance to the results obtained* 

6*2 Land Bse Map t 

*Slm land-us®' map prepared in the paeogra® BCL^SS 
is in fact a mosaic prepared from separate ccwpputer print 
outs* Each coBspjttfflc printout consists c£ aa overstrucJ©- 
character representation of C 120x79) pixels* Six such 


I: 
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sheets were pasted to get a laosaJxs and four, sucli aweiaics 
covered whole of the study area# The relative positlQBis 
of iiKJSaics, id<»itified as A,B,C and^ D, are represented In 
Fig i6*x)m h moealc wise discussion of them is given in 
one o£ the following secticxis* 

6»3 Field Suarvev * 

To check the accuracy of prepared land »se map an 
esKtensive field siavey was conducted. *Grc«ind truths* 
were collected during this survey# This survey was carried 
out in month of WiarCh* The data collected during this 
period is given in table form# This <3feta w^ craeapared 
with the one furnished 1^ the assumed land use map 
a ccMitingesaKjy table was calculated to give the accuracy 
in percentage for different classes.# The table is given 
in Table {5#2)# 

e*4 Discussion over the tend Ose Ma.n *» Mooealc k s 

The »ME»st strickingly identifiable pattern on this 
»osaic Is the meander of river Gomati# Most of the area 

covered is cultivated# with study patches of builtrp areas# 

\ 

The course of river has not Changed much since 1972# the 
year in ^Ich the Sirrvey of .India released ^the topoebeet 
of the' area# This statement is sus^jorted the prepared' 
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land use map as ymtl as the ^grcusid truth * ■sunieisf*. saiicl 
bar fra» 3i.ite nwSimc 1859 pixel liiasber 1570 to line ISta3(ft^sQr 
1864 pixsl hixiaer 1545 Is a recjetit featisfe hot (depleted In 
top<^heets« It se^wt that dariag high discharge ^s^soh: a 
subsidiary chaohel trapesed the path -and dried aftsar r« 2 e« 
ssioEi of floods leaving the existing sand bar* ^t' siost of 
the places the builti^ area was correctly niapped (4 out 
5 sites). The laiidi use laap shonrs patches of sand in betwe^ 
agricultural land* Certainly^ it is very uhlihely that saad 
be there, hence few such sites w€3re visited Airing the field 
aurvciy* It was f ot»d that these parts of land ba«s not been 
used fesr agricultural purposes for seny years because of its 
piDor productivity* Siot tilling the laynd for raany years co®* 
pacted ttte soil and tl^ ccsi^^cted soil, which is st^jposed to 
has the sa^ spectral respmse a&.deryBe sand was classified 
as sand* Thotsgh this class of unproductive land was not 
taJeen into consideration initially »hi2® selecting training 
data set, it was a««»pted as a sespasrate class after the 
field survey* The deeper vatesr of river ^Sanga is rightly 
classified and. shown in iRictel.l« of ri^t hand side of the 
mosaic* The road C5JH29> and i^iljway Itoe RallwaysCJhapaaj 
Varanasi, Allahabad Br^sdi) could -iiiot be classified because' ci 



coarse resolution of lia^lsat* Ho latentif laifele growtlj of . 
built up area toolfe place In ecoCemed area# 

Mosaic B » 

ihls TOsaic covi*« parts, of Varanasi aii€ Obasiptir 
districts# *lhe curvilinear patteibEi tbe river Oanga 
is very truthfully reproduced in land use nap in accordance 
with what is there in the toposfeeet ( SOl, 63 .0/3 ) eaccflqgrt: 
slight changes in the shape of the sand islands in river 
water# Ihe sire of sand inland (line nuBP*>er 1931 pixel 
nuBtoer 1670) has grown bigger as coai»red to cne given 
in the toposheet due to deposition# Siiailarly the/ shape 
of the sand islaai^ (line ninter 1970 pixel nunher 1746) 
had also change# The Imiltii^ a^a shewn here m& there 
on the islands are the constructed banks (ghats >• 

the change in builtt^ emm is vesy premlnesit in 
this nosaic when ccMi^pared to the status in 1972# Two 
towM have shown noticeable growth in Isailtiip artsi '•* 

(1) abivdaspur Gtoagabarar (line niwb®r 1950i pixel mmheac ^ 
i66S) and (2) ^hnrn Har^wa (line nunbeer %&snr pixel number 
1705)# Chanda Tal which is asidy pend is al»0 Clearly 
classified as less doinh water class Clii^ nunfeer 1^0 
pixel nuBtoer 1709)# Soieeti«»s roaite and railway lines, 
Inspite of their tfelckaess less^itlian the s!S»tial resolutlc® 
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o£ MSS, appear very qlearly. It is because of tbe builtup 
area coraixig along tbe roads and railway lines* Hiey sfaow 
a distinctive linear pattern of builtup areas as is evident 
from Une nunftjer 1803 pixel ntaaber 1756 to line nisitjer 1879 
pixel nimdaer 1826. 

Uiere is a strilcing dbange in cc«if luenee point o£ 
river Ccxaati and river Oar^* As tbe toixjsfltieet 63 j shows 
that the confluence of these two rivers was at Saldpur (line 
number 1790 pixel nui&Msr 1691) In 1972* But the lm& use 
map that was prepared Showed a deviation of confluence to 
a point upstream Cline number 1829 piNrei nuatoer 1670)* the 
results of field survey will W present later on* 

Mgpaic C • 

this mosaic c«xtsists of land use map covering most 
of Varaawasi city* Again the riv«t Oanga is very prominentely 
and accurately mapped* the meander patterns remained same 
at least for last 16 years*, only the ^pe and siae of saaad 
hfsivs have chaoiged* "to the leest of river <3anga lies the city 
^ose northern limit is the riv®r v&xtxtmp the accuracy with 
which the city is mapped is more than agreeable* Xt was 
chedked field Bwcv&y* 'Hhe conflissnce of river TanHia 

with river Ganga could not be mapped probably owing to its 
smaU thicteess* 
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On the east oE riLveir Oanga lies saall toims as ■ 
Bahadurpur CliJic nuiaber 2097 pixel niiatjer 1553) , Saia* 
nagar (line nunOser 2155 pixel fti]ni>er 1540), 
colony (pixel niMser 1599 line nuiatoer 2135) and villages 
as Scsrnra (line nOTt>er 2128 pixel nnitiber 1516), Ocrari 
(Une nwstoer 2115 pixel mraber 1522) etc* The rasp also 
contains Grand 'Eruiik road and m Railway loain line 
because they appear ccmtlnuously in raany pixels* The 
OT road is the upper dlagcxial line and 13®, naain line 
Is the Icwer diagonal line in the right half osE the rasp* 
Pcxid near Wajldpur (line nuraSber 2149 pixel iiiraasar 1545) 
is also rightly classified as water oE less depth* 

•Ihe boom oE development of bulltup area is felt 
in locality covering towns Kunda Kalsn., Bhawanipur, 

Narla (.line ntratjer 2093 pixel msd&er 1600)* It was 
found to be true during the field surrey* 

l4o«aie D * » 

This raoaaie covers the town of fitaghal Saral* In 
the upper part c£ inosaic b«®d of river GSsuaga is shown 
along with a low discharge dtmuam 1 separating from the 
same hnown as Sota nala#' which. ta3«es a turn in S shape 
to li^t ^e pain branch e£ Oanga rivssr* The bulltup 



area near Jalupur (line nuaiber 2004 plrel nuab^ 1677) 
showed an increase# In ftighal Sar8|.i NR main ?-| , ne a ieet a 
the NR main line {Electrified)# This clu^turiUag c£. 
railway lii:^ produces a misls^-ding classiflcatioii# The 
reflectance walues p-irely from a railway track se«@ to 
be lower wh€«i ccw^rfKl to the s-sro tsS. built up area m: 
qultiirai^ land# ^ere there is a single track the 
spectral respcaise# reccarded for a pixel cowering that 
area# is gowemed laostly by mp&ctsml rehouse of eulti* 
wated land and/or built up area close to railway track 
and that pixel is classified as either built \sp area or 
cultivated land# But the responses, for pixels representing 
multitrack patch of land# the r^ponses of raihmy line 
predominate# Since the lacightness values of .railway track 
are low as ccsmpared to that of built up area# such pixels 
representing wilti<»track railway lines are eiassliied 
into a c lass that has a lower l^ight^ss an as compared 
to that of ^l>uilt up area i© water of less deptli^ Ibis 
®d.cclasol£lcatioo~«cises'out’a£ tiie fact lhat less nu»i>er 
of ciacc "‘s# to which a pixel is to be assigned# are consi- 
dered becMUSe follcwi^ two facts - ,<l)..the resolutic® 

Is coarse hence (MJ@ can not always expect to discri- 
minate all the classes present 'em the ground ^ (2) im- 
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avall^lllty of detailed (large scale) roaps for sel« 2 tiaEi 
of training data points fear all the classes^ Siailar 
phenomenon of isaisolasslficatioii oectired in mosaic C also* 
Orand TTimk Road (I3H2) is very clearly mapped in the mosaic 
D# Bnilt ijp areas have started coming np in sarro«®ding8 
of Shivadasa Gangbarar (pheel nimdaer 167©, line nnstier 1945 >• 

6*5 Acenraev Asaesment t 

Based cm field survey contingency table (T^le 6*2) 
was prepared f or accuracy assesment. Both the errors - 
errors of omission and error of camnission, are calculated* 
The table is self eaq^lainajtory* The accuracy is glv®a as 
percentage* 

6*6 CtxKslusion I 

The automated digital approach to 1|SS| use map 
preparation is a worth while process, in s^se that the 
whole process starting from the data acquisiti^raa, data 
correctic® and information «rtractic«i finises in a short 
time spip* Where as the traditicnal methods of ground 
surveys and map preparation tid»s as mixsh as ...two years t j 
Another benefit that can be incurred from the me of 
satellite plateforms is tl^ fast and repetSvf data about 
the land cover and lai»3 use* The use of i©s landsat 



data with a coarse resoluticoi of 80 ra have given a tolerable 
accuracy, few of the level II classes were also effeciently 
classified* Fowc land use classes namely water, compacted 
sand, asgricultural land, and built up area with fto-ther 
c lassif icatic^ oi class water in two level II classes has 
been achieved* 

6*7 Further ReecBtamendatioiis t 

The resoluticm of lendsat MSS is about 80 at* Wid» 
such coarse resolution the classification at level II is 
not very fruitful* Thus good resoluticm data as from 
Landsat TM or SPOT should be tased for effective and 
reliable level II classification* 

The spatial change in the laijd use pattern over a 
span of time can be studied, if landsat data are available 
at two tijoes* Ibis stu% should be ecxidiKJt«l in future 
toy scaaeone interest®3| by aeguiring a new landsat imag'i&rf 
of same area* ' 
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Ground Truth Siarvey results* 




No* visits 




use lama, use 
iccaticn cJbss on olass on 
site ma' 


BeE^rks 


h Malahia 


Madarwa 


Bhogpura 


Btielupura 


E Sarian 


F Delwaria 


25*14 *50”N 
83* 1*25*E 


25*15 *16 "N 
83* 1*10*E 


25*17 « 0"N 
83* 0*10**B 


25*17 ‘SO^N 
83* 0*20*E 


25*20» 3”N 
83* 1*27''E 


25*20»25"N 
83* 0*10*E 


Correct classic 
Correct classlf 


Coexist classlf 
Nrc*]g classif 
Correct classif 


Correct classic 
Wrong classif 


Correct classif 
CcBcxect classif 
Correct classif 


Correct classif 
Ccorect classif 
Pftroog classif 
Ccacrect classif 


Correct classif 
WrcKig classif 
■^eaig classif 


6 Seesari 


25*17 ‘SO-N 
83* 2*i7*S 


Corr®st classif 
Correct classif 
Correct classif 
Cccreet classif 


(contd*** ) 



H JalilpixT 

25*18 *15«N 





83* 3* E 




22 


1 

1 

Correct cdasslf 

23 


1 

52 

Wrcaig classlf 

24 


2 

2 

Ccarrect classlf 

25 


52 

52 

Correct classlf 

1 Raiaiagar 

25*16 »20''N 





83* 2* 5”E 




26 


1 

1 

Gcrri^t classlf 

27 


1 

1 


28 


2 

2 

m 

29 


2 

2 

m 

30 


52 

52 

1 

J Sholapur 

25*16 • N 





83* 4*10*1 




31 


1 

1 

■tf' 

32 


2 

2 

If 

33 


2 

52 

Wrong classlf 

34 


77 

77 

Correct c lasslf 

35 


2 

77 

Wrong classlf 

K aaghalSaral25 *15 *10"R 





83* 5 *12“! 




36 


1 

1 

Correct classlf 

37 


2 

1 

Wrong classlf 1 

38 


2 

1 

Correct '’claOSlf 

39 


1 

52 

Wrorg classlf [ 

40 


52 

52 

Correct classlf i 

41 


1 

77 

Wrong classlf 

42 


77 

77 

Ccsrreot classlf ■ 

43 


1 

t 

Wrong classlf 

44 


77 

77 

' Correct classlf 

h Sallnabad 

25*19* 5"ll 





83*10*10*1 




4S 


1 

2 

Wrong classlf 



2 

2 

Correct classlf 

47 


1 

1 

CcRcrect claussif , 





CcO£ltd«*«) ! 



5 


M Ifeyrilcot 

25*16* 3*W 





83*14 * E 




48 


1 

1 

Correct C-lassif 

49 


2 

2 

Ccsrrect classif 

50 


52 

2 

Wrong classif 

N Gokulpura 

25*19 •40'*N 





83* 9»35"E 




51 


51 

51 

Correct classif 

52 


52 

52 

m 

53 


2 

2 

m 

54 


73 

73 

m 

56 


52 

52 

m 

0 Raynagaxha 

25*25 'ICN 





83* S’^O^E 




57 


1 

1 

Correct classif 

58 


2 

2 

ff 

59 


rt 

2 

Wrong classif 

60 


77 

77 

Correct classif 

61 


52 

2 

Wrong classif 

E Sarnath 

25*18 •10'*N 





83* 2* 




62 


1 

1 

Correct cl»sif 

63 


1 

1 

n 

64 


2 

2 

'll 

65 


2 

77 

Wrong classif 

66 


Ti 

77 

Correct classif 
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67 

68 

69 

70 

71 
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Khtjard 


25*27*45"E 
83* 2*30''E 


1 

2 

2 

77 

77 


2 Wrong classif 

,2 Cocreefe ciassU 

17 Wrong classif 

2 WtcMfig classif 

'77 Correct classif 


Correct classif 

m 

m 


72 

73 

74 


2 

2 

2 


2 

2 

2 


IconttS* • • ) 


75 


2 

2 

Ccarrect 

76 
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■il' 

77 
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77 
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79 
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2 
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80 


77 

77 
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52 

52 

m 
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52 

52 

^ ■ 

83 


1 

1 

Correct 

84 


1 

1 

•11 

35 


52 

52 

tl 

86 


77 

Tt 

n 

87 


2 

77 
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88 


1 

2 


89 
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77 

2 

tl 

90 


1 

1 
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91 


51 

51 

if 

92 


52 

52 


93 


52 

52 

n 

94 

V *«iawari 

25*24* 5*11 
83*11*10*8 

1 

2 

wros© 

95 


1 

1 

Correct 

96 


2 

2 

m ■ . , 

97 


51 

51 

^ m- 

98 


73 

73 

■:S' ^ 

99 


73 

2 

WlETOlig 

lOO 


52 

52 

Correct 

101 

w Jtlsirpw 

25*24 *4S*H 
83* 9*40*K 

1 

1 

IRt 

'f fit ' 


' 1 ■ 

^ t 

Wrong cla®el£ 


(cC 81 t<l* # » » ) 




103 

104 

105 

106 

X Ajgara 25*28*30"H 

83*12*15“! 

107 

108 

109 

110 
111 
112 

113 

Y Tanda 25*30 * 0"N 

Kalan 83*10* 5*1 

114 

115 

116 

117 

118 

119 

120 

2 Saldpura 25*32* 5“W 
83*16* 0“! 


4 5 6 


52 

52 

Correet 

2 

2 

m 

51 

51 

li 

73 

73 

II 


1 

1 

Ccarrect 

1 

2 

it 

52 

2 

Wrong 

2 

2 

Correct 

2 

77 

I'lTotig 

77 

77 

Correct 

1 

1 

m ' 
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1 


73 

73 

#1 

51 

51 

n 

2 

2 

m 

2 

% 

m 

52 

52 

m 

73 

73 

m 


121 

122 

123 

124 

125 
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127 


t 

1 

m 

51 

51 

m 

51 

51 

m 

52 

52 

II 

2 

Tt 

Wrong 

77 

77 

Correct 

2 

77 

Wrong 


78 


o m m 

m 


tHf 



s 

'* 

♦ 

• 


o 

CM 

m 

<M 


CM 

r4 

01 

m 

m 


f*H| 



W 

I 

f 

f 

m 


w i O 


4J If & 

g 

^ l<-% S«l«4 

HI gN •»» »tA 

9 J> 
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APPENDIX TO CHAPTER 6 


This appendix ccmtalns the coosale prepaired fcxr the 
study area and results of Image paroc^slxig of the dlglti* 
zed data for the same* 


Plate 1 I shows the image of the Band 1 data* Since the 
Brlglptness values were reqcarded frc»a 6 to 130# hence does 
not use the full rai^e available cai IP syst^# ie O to 255* 
It is the reascw for Its pocar contrast result* 


Plate 2 * iftiows a stretched image of a part of the sceaae 
in band 5* Hie ^leari^ isf the image shews the advantage 
of image stretching over unstrelKihed one (as in plate 1 )• 


Plate 3 # is a result of superposition of stretched data 
in three different bae^s# ie band 4# band S ai^ band 7* 

It is zoomed to bring out certain features which were not 
iCNsitified on the llne»*prlnter map because of pjcr reso-* 
lutionr but the spectral resoluticm as well as the spatial 
resolution of IP system being good# it brwaght out many 
features prominently as the roads# railways et®* 

Plate 4 t it is again a^ zoomed# superimposed .and stretched 
image In three channels* Blt^ colom represents deep waters 
as in river alixig its c^tre line and few ponds* Grayish 



so 


r«a »how« the cultivattsd area 'when a$ red dtiaiitiel represeaats 
sand bars and unused barren lardir Black colciJir represents 
water o£ 1mm depth as well as the built mm&t because 
it is not possible with eyes alcxie to differentiate between 
slightly varying tcme these two classesn Ihe green colour 
represents the built area made mit asid thatch that 

have less reflectance as conapai^d to the ccncrete or raasjtonary 
structure# 


» 
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